Quantum behavior of FRW radiation-filled universes 
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We study the quantum vacuum fluctuations around closed Friedmann-Robertson-Walker (FRW) 
radiation-filled universes with nonvanishing cosmological constant. These vacuum fluctuations are 
represented by a conformally coupled massive scalar field and are treated in the lowest order of per- 
turbation theory. In the semiclassical approximation, the perturbations are governed by differential 
equations which, properly linearized, become generalized Lame equations. The wave function thus 
obtained must satisfy appropriate regularity conditions which ensure its finiteness for every field 
configuration. We apply these results to asymptotically anti de-Sitter Euclidean wormhole space- 
times and show that there is no catastrophic particle creation in the Euclidean region, which would 
lead to divergences of the wave function. 
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I. INTRODUCTION 

Homogeneity and isotropy of the universe on large scale 
is a good approximation to describe the classical behav- 
ior of the universe. Friedmann-Robertson-Walker (FRW) 
models are specially designed to implement these proper- 
ties. Nevertheless, seeds of inhomogeneity and anisotropy 
are needed in order to describe the cosmic structure. For 
this purpose, studies of cosmological perturbations are 
necessary. Seminal works in this direction were done in 
Ref. [1] and later on in Ref. [2] where the authors studied 
the stability of de Sitter space. 

The classical description of the universe breaks down 
for energies of the order or above the Planck scale. There- 
fore, it is necessary to use a quantum theory of gravity 
and to postulate some boundary conditions for the uni- 
verse in order to describe its initial state. Despite the 
absence of a fully consistent quantum theory of gravity, 
many studies have been carried out that shed light on 
the problem of the creation of the universe with differ- 
ent boundary conditions [3-8] . These works characterize 
the quantum behavior of the universe in the semiclassical 
approximation through its wave function in both min- 
isuperspace and superspace, where the inhomogeneous 
and anisotropic modes are included perturbatively in the 
models. 

In Ref. [9], it was noted that the wave function of a 
closed FRW universe with a positive cosmological con- 
stant becomes infinite in the forbidden (tunneling) region 
when the universe is filled with radiation and subject to 
vacuum fluctuations of a massive scalar field conformally 
coupled to gravity. In other words, the author concluded 
that during the tunneling process a catastrophic parti- 
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cle creation takes place. He also speculated that perhaps 
these phenomena might be a rather common feature of 
tunneling processes due to quantum gravity effects. In 
Ref. [7] , it was shown that the wave function of a de Sit- 
ter universe in the presence of gravitational perturbations 
increase for some boundary conditions but never diverge. 
Similar results were obtained in Ref. [8] for a minimally 
coupled scalar field and tunneling boundary conditions. 

In this paper, we develop a method based on Ref. [8] 
to study the quantum behavior of the wave function of 
a radiation-filled FRW universe with cosmological con- 
stant and radiation, which includes vacuum fluctuations 
represented by a massive scalar field conformally coupled 
to gravity. These vacuum fluctuations will be regarded 
as perturbations to the homogeneous and isotropic solu- 
tions of the Wheeler-DeWitt equation. We can deduce, at 
least for some values of the scalar field mass and negative 
cosmological constant, that the perturbed wave function 
is not divergent in the classically forbidden regime. As 
we will see, the finiteness of the wave function is due to 
the regularity and boundary conditions, which although 
restrictive still allow for finite solutions. These quantum 
states represent asymptotically anti de Sitter wormholes 
[10,11]. 

The paper is organized as follows. In Sec. II, we re- 
view the classical behavior of a closed radiation-filled 
FRW universe with a cosmological constant, both in the 
Lorentzian and Euclidean regions. In Sec. Ill, we derive 
the Wheeler-Dewitt equation for these universes in the 
presence of vacuum fluctuations of a conformally cou- 
pled massive scalar field and perform the semiclassical 
approximation. In Sec. IV, we deduce the general match- 
ing conditions that relate the wave function defined in the 
different semiclassical regions. We also impose the reg- 
ularity conditions. In Sec. V, we obtain the background 
wave function and linearize the equations for the mat- 
ter vacuum fluctuations thus obtaining generalized Lame 
equations. We solve these equations for asymptotically 
anti-de Sitter wormhole spacetimes. We show that this 
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perturbed wave function is finite for all possible values of 
the scale factor and scalar field configurations. Finally in 
Sec. VI we summarize our results and conclude. 



II. LORENTZIAN AND EUCLIDEAN 
BEHAVIOR OF FRW UNIVERSES 

The main part of this paper will be devoted to study 
the quantum behavior of a closed FRW universe filled 
with radiation against perturbations due to a massive 
scalar field conformally coupled to gravity. But before 
that, let us shortly review the classical behavior of a 
closed FRW universe filled with radiation [12]. In our 
analysis we include a cosmological constant A = 3A, and 
we represent, for simplicity, the radiation of our universe 
by a conformal scalar field A. The FRW metric can be 
written as 

ds 2 = a(fj) 2 (-dfj 2 + dSlj), 

where fj is the Lorentzian conformal time and dft 2 is the 
line element on the unit three-sphere. Writing the radi- 
ation field as 



A(x,r?) = 



the Lorentzian equation of motion for x{v) becomes 

X = 0, (2.1) 



x" 



where the prime denotes derivative with respect to fj. 
The equation for this field can be integrated to obtain a 
constant of motion K, related to the energy density by 
p = K/a 4 : 



x' 2 + x 2 



3K 
2G' 



where G is the gravitational constant. Then the scale 
factor must satisfy the equation 



'2 



+ V(a)-K = 0, 



where 



V(o) - a 2 - Xa 4 



(2.2) 



(2.3) 



The shape of the potential V(a) depends on the sign of 
the cosmological constant A. For a positive cosmological 
constant, it increases up to a maximum value at a = 
^= , and decreases after that for scale factors larger than 
-^7= Fig. (1). For a negative cosmological constant the 
situation is rather different, as the potential is always 
increasing and never negative Fig. (1). 

We can distinguish three kinds of behavior for a. The 
first one describes a collapsing universe. This is the case 
when the cosmological constant is negative and K ^ 0, 
for which 



a{rj) = a p cn 



1 



2m 



(2.4) 
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FIG. 1: This figure shows the potencial V(a) defined on 
Eq. (2.3). The darkest and lightest curves correspond to a 
positive and negative cosmological constant A = 3A cases, 
respectively. The parameter K is related to the amount of 
radiation presents in the FRW universe. On the one hand, 
for positive A and K smaller than the maximum of V(a), as 
in the case plotted, a_ represents the maximum radius of the 
collapsing radiation-filled FRW universe, while a + represents 
the minimum scale factor of the asymptotically de Sitter uni- 
verse. On the other hand, for negative A the scale factor a p 
represents the maximum radius of the collapsing universe. 



where rj <G [r)-,r) + ], a 2 p 
and y/m(f)+ - fj-) = K( 2 , m 



= m = y/l - 4KX 
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-). In these expressions, 



cn[x, is a Jacobian Elliptic function and K( I ^L) i s 
the complete Jacobian elliptic integral or quarter-period 
function [16, 17]. Note that fj + is an arbitrary constant 
that can be set equal to zero. The scale factor of this 
universe increase from a = at fj = fj- up to a = a p 
for fj — 0, which is the maximum radius of this universe. 
The other case of a collapsing closed FRW universe cor- 
responds to a positive cosmological constant and a value 
of the parameter K, related to the amount of radiation 
present in the universe, smaller than the maximum of 
the potential V(a), i.e., < K < j^. Under these condi- 
tions, the scale factor in terms of the cosmological time, 
dt = a(fj)dfj, has the expression 
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2A 



a(t) 2 = — < 1 — m cosh 



2\/A(i-i_) | 



the maximum value of the scale factor is a 2 _ = , the 
corresponding to i = t- which is the solution of the alge- 
braic equation V(a) = K. For both solutions, the maxi- 
mum radius of the universe increases with the amount of 
radiation, given by K. 

The second kind of solutions describes an asymptoti- 
cally de Sitter space-time when A > 0, whose scale factor 
is given by 

a(t) 2 = ^-{V 4:KX sinh \2V\(i - i+) - 
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— cosh 



2>/X(t-t+)] + l} 



for K > and 
a(t) 2 



1 

2A 



1 + m cosh 



2>/A(t - t+)] } 



for if < jj, where t € [i + , +00). The difference between 
these two cases is that for sufficient radiation, i.e. K 
is larger than the maximum of the potential V(a), the 
scale factor grows from zero up to infinity, while in the 
opposite case; i.e. for K smaller than the maximum of the 
potential V(a), the scale factor grows from a minimum 
value different from zero, = , due to the presence 
of the barrier of potential V(a), to become asymptotically 
de Sitter. 

And finally, there is a third kind of solutions which 
exactly coincides with a de Sitter space-time 



a(t) 2 — i cosh 2 
A 



V\(t-t+) 



in the absence of radiation and for a positive A. 

It can be checked that there are not classical solutions 
of the Einstein equations corresponding to a closed FRW 
in absence of radiation, K = 0, with a negative cosmo- 
logical constant. It is only for K > that it is possible 
to have a Lorentzian evolution for the scale factor a. 

Up to now, we have described the different possi- 
ble Lorentzian solutions for a closed homogeneous and 
isotropic universe filled with radiation and we have seen 
that the potential V(a) forbids the classical evolution for 
some values of the scale factor. Therefore two classical 
FRW universes filled with radiation with K < -A- are dis- 

4A 

connected and the scale factor for a FRW universe with 
a negative cosmological constant has a maximum value 
when the content of the universe corresponds to radia- 
tion. 

As it is well known, the fact that two classically al- 
lowed universes separated by a potential barrier are clas- 
sically disconnected does not mean that they cannot be 
connected quantum mechanically. In the lowest approx- 
imation, this connection is established by an instanton 
whose explicit form can be obtained by performing an 
analytical continuation of Eq. (2.2), for a positive A, so 
that the classically forbidden region is now the permitted 
one. The solution for the scale factor must satisfy 



a{r}-) = a_ 



a{v+) = a+, 



in order to connect with the two classical FRW universes. 
From the analytically continued version of Eq. (2.2), we 
obtain the following solution for the scale factor: 



1 \2 1 + ™ , 2 

<V) 2 = dn 2 



1 



'■(v-v+), 



2m 



m 



(2.5) 



where 77 e [V-,V+] with 77+ - r?_ = -y^^K( I 2 ^). In 

this expression dn[x, jzp^] is the Jacobian elliptic delta- 
amplitude function [16, 17]. Note that 77+ is an arbitrary 



constant that can be set equal to zero. This instanton 
was also found in Rcf. [12], where the authors considered 
a closed FRW with a material content corresponding to 
a massless scalar field conformally coupled to gravity. In 
the absence of radiation, m = 1, the turning points of 
the potential V(a), i.e., the solution of V(a) = K, (sec 
Eq. (2.3)) becomes a_ = 0, a + = \/T/X. The instanton 
(2.5) then acquires the simple form 

a (v) 2 = jSCch 2 (T]-T] + ) 

and 77 + — ?7_ = 00. 

While for a positive cosmological constant the Eu- 
clidean solution for a closed FRW universe filled with 
radiation connects two classical solutions, for a negative 
A the solution behaves as an Euclidean asymptotically 
anti de Sitter wormhole. This can be easily deduced 
from the analytical continuation of Eq. (2.2) to imagi- 
nary conformal time. The scale factor of the wormhole 
looks like [10] 



a(r]) 2 = o 2 nc 2 



Vm(f] - 77-) 



1 



in 



2m 



(2.6) 



where 77 G [??_, 77 + ] and ^/m(r] + — t/_) = K(^i ? ii). In this 
expression nc[x, 3^] is a Jacobian Elliptic function [16, 
17]. The value a p describes the radius of the wormhole 
throat. 



III. THE WAVE FUNCTION OF THE 
UNIVERSE 

The quantum behavior of the FRW universe can be de- 
scribed by the solution to the Whceler-DeWitt equation 
[13]. In the WBK approximation, the wave functions 
can be approximated, under certain conditions, by in- 
going and outgoing modes defined through the classical 
action in the Lorentzian section, while in the Euclidean 
sector, the wave function can be approximated by lin- 
ear combinations of increasing and decreasing modes in 
terms of the Euclidean action. Boundary conditions that 
determine these linear combinations are also necessary. 
In this section, we will obtain the general shape of the 
wave function of a closed FRW universe filled with radi- 
ation and whose content corresponds to a massive scalar 
field conformally coupled to gravity. 



A. Canonical formulation 

We will consider a minisuperspace described by two 
degrees of freedom, the scale factor a and a homogeneous 
and isotropic scalar field conformally coupled to gravity 
A. Around this minisuperspace, we will study the linear 
perturbations due to an inhomogeneous and anisotropic 
massive scalar field $ conformally coupled to gravity. We 
will obtain the Wheeler-DeWitt equation from a specific 
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representation of the Hamiltonian of the system which 
can be constructed easily from the classical action of the 
system 

S = >S g -+- S Y -\- S m , 



l — ( d 3 xVhK, 



s r = 



+ 



8ttG 

d 4 x 



-(VA) 2 



12 



-i?A 2 



+ 



i J d 3 xVhKA 2 



d 4 x. 



+ 



i J d 3 xVhK<S> 2 , 



+ 



(3.1) 



where G is the gravitational constant, /i is the mass of the 
scalar field and K is the trace of the extrinsic curvature. 
We have used the sign of conventions of Misner, Thorne 
and Wheeler [14]. 

We introduce new variables which correspond to an ex- 
pansion in hypersphcrical harmonics of the massive scalar 
field around the background solution described in the 
previous section as follows 



$(x,77) = a(r?) 1 



<P(v) + J2f^)Q nlm (x) (3-2) 

nlm 

m)+Y.^{fi)Q nlm ^) 

(3.3) 



nlm 



where Q nlm are the scalar hyperspherical harmonics, 
eigenfunctions of the 3-dimensional Laplacian V 2 in 
the three-sphere, i.e., they satisfy the eigenvalue equa- 
tion \7 2 Q nlm = -n(n + 2)Q nlm , with n = 0,1,2,.... 
The mode n = corresponds to the homogeneous and 
isotropic perturbation, while higher values of n corre- 
spond to inhomogeneous and anisotropic modes. We 
have considered the background solution of the massive 
scalar field $, equal to zero, i.e. if — 0. From now on and 
for the sake of simplicity we will drop the indices Im and 
keep only the eigenvalue index n in all the expressions. 

The space-time metric must also be expanded around 
the homogeneous and isotropic background solution. If 
we write the Lorentzian metric in the 3 + 1 form 

ds 2 = a 2 (-N 2 + N'N^drj 2 + 2aN l dridx l + g l0 dx l dx\ 

this expansion can be written as 

N = 1 + ^T^Q", 

n 



(a n Q? 3 + b n PP 3 + c n S?j + d n Gl) 



Slij being the metric in the unit three-sphere and Q n , P n , 
S n and G n are the standard hypersphcrical harmonics on 
the three-sphere [1, 6]. 

The coefficients x n = {a n ,b n ,c n , g n , j n ,k n } can in 
principle be eliminated by means of a diffeomorphism 
on the three-sphere and choosing suitable lapse and shift 
functions [6, 7]. The only terms in these expansions that 
cannot be gauged away correspond to pure transverse 
traceless tensor perturbations that describe gravitational 
waves. These are represented by the coefficients d n . The 
Lorentzian action up to second order in perturbations 
has the form 

S = S {0) [a,x\ + sfl r [a,x,Cn-,x n ,d n ] + Sf£>[a,f n ]. 

This can be seen from Eq. (3.1) taking into account that 
the massive background field vanishes; i.e. ip = 0. In- 
deed in this case the perturbations of the massive scalar 
field decouple from perturbations of the metric and the 
radiation up to second order. Since we are interested in 
the behavior of closed FRW universes against the quan- 
tum fluctuations of the vacuum of the massive scalar field 
$, we have chosen its homogeneous background mode to 
vanish. So the explicit expression of S g + r is not neces- 
sary to study the behavior of the FRW universe under 
perturbations of the massive scalar field. The zero order 
action and the second order action of the scalar field 

(2) 

perturbations S m have the form 



3tt 
4G 



( a 2 -a' 2 -Aa 4 ) + ^(x' 2 -X 2 ) 



Jr)- 



(3.4) 



where 



U n (fi) = (n + l) 2 + f i 2 a(fj) 2 



Finally the Hamiltonian of the system can be written as 
H = H^[a,x] + H^[a,f n }, 



(3.5) 



This Hamiltonian describes the classical constraint of 
our system H — and is related with the invariance 
of the Lorentzian action under time reparametrizations. 
This constraint will be our starting point for studying the 
quantum behavior of closed FRW universes filled with 
radiation. 

The constraint H = in the context of quantum grav- 
ity becomes a constraint on the wave function of the uni- 
verse leading to the Wheeler-DeWitt equation, which can 



■5 



be written as 

G_ 

3^ 



^;+^(« 2 -Aa 4 ) + ^(^-x 2 ) + 



perturbative part of the wave function of the universe. 
Using the WBK approximation, we obtain 



1 + °° 1 ~ 



where 



£/„(a) = (n+l) 2 + M V. 



The functional dependence of the wave function on the 
radiation field can be obtained by separation of variables. 
The part of the wave function *(a, /„) which depends on 
the other degrees of freedom present in our model, a and 
/„, must satisfy 



■^ + ^Vo(a) + Y,(-dl-V n (a).f n 



71=1 



*(a,/„) = 0, 
(3.7) 



where the potentials V n (a) and V (a) are defined as fol- 
lows: 



V n (a) = -[(n + l) 2 + M V], n= 1,2,. 



y (a) 



a 2 - Aa 4 - X. 



(3.8) 



Here K is a separation constant, related to the energy 
of the mode x an d quantifies the amount of radiation 
present in the universe as we have seen in Sec. II. 



fdSoY _ 9tt_ : 
\da J ~ 4 
dSp dS n 37r 2 _ 
~da~~da 2 n ~ 



V (a), 



-3ttK, 



(3.10) 
(3.11) 



where we have performed an asymptotic expansion on G 
and likewise for Rq and R n . Similarly to the usual WBK 
this approximation is valid as long as 



G 



dV (a) 
da 1 



< |Vb(a)|*. 



(3.12) 



From now on, we will consider that the amount of radia- 
tion present on the closed FRW universes in the case of 
positive cosmological constant A, is such that this kind 
of models has two classically disconnected solutions, i.e. 
< K < 1/4A, allowing for quantum tunclling between 
the two universes. As can be seen from the validity of the 
WBK approximation Eq. (3.12), this condition breaks 
down near the points a_ and a + in the case of A > 0. 
These points corresponds respectively to the maximum 
radius of the collapsing universe and the minimum radius 
of the asymptotically de Sitter universe. Similarly the 
WBK method is not valid near the point a p which corre- 
sponds to the maximum radius of the collapsing universe 
in the case of a negative A. Therefore we have to ana- 
lyze the behavior of the wave function near these turning 
points using other methods as explained below. 



B. Semiclassical approximation 



C. Turning points 



We will treat the quantum behavior of radiation-filled 
FRW in the semiclassical approximation [7, 8] , where the 
physical lengths involved in our problem arc larger than 
Planck length l p . In this approximation, the solutions to 
equation (3.7) will be written as linear combinations of 



*i(a,/n) = cxp 

*2(0,/n) = CXp 



I 1 + °° 

n— 1 



(3.9) 



where and ^ will be decreasing and increasing func- 
tions in the classically forbidden regime or ingoing and 
outgoing waves in the classically allowed regime [8]. Wc 
will deal firstly with the regions of validity of the semi- 
classical approximation in term of the amount of radi- 
ation present in the universe and in later stage of our 
study, we will specify clearly the kind of functions on 
each regime, they will depend on the inclusion of the 
Lorentzian or the Euclidean time. So and Rq are related 
to the unperturbed part of the wave function of the uni- 
verse while S n and R n are related to the vacuum fluctu- 
ations of the massive scalar field and corresponds to the 



Starting from the expansion of the potentials Vo(a) 
and V n (a) around each turning point, a_, a + and a p , the 
Wheeler-DeWitt equation (3.7) acquires the form 



G 



37r dVo(a) 



3tt a 4G da 



\a= ai {a-ai) 



+ 00 

£ 

71=1 



1 



-d) n V n ( ai )f n 



*(a,/ n )=0, (3.13) 



where a, = a±,a p . This linear approximation hold if and 
only if: 



d 2 a 



da 



(3.14) 



The wave functions of the closed FRW universes can 
be expressed as linear combinations of \I/i(a, /„) and 
^2{a, f n ) defined in Eq. (3.9), where So, Rq, S n and 
R n satisfy the differential equations (3.10,3.11), if the 
value of the scale factor, a, is such that the condition 
(3.12) holds. While in the linear regime the behavior of 
the wave functions of these universes will be solution of 
Eq. (3.13). As we will see there are values of the radius of 
the universe where both conditions Eqs. (3.12,3.14) hold 
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and it is possible to connect the wave function on the 
WBK regimes through the linear regime. 

Let us begin analyzing the case of a positive cosmo- 
logical constant, in which the conditions for the linear 
approximation are 



o+| « 



2m 



2 + 3m 



a+, \a- 



2m 



-2 + 3m 1 



(3.15) 



respectively around a+ and a_ , while near these turning 
points the WBK conditions (3.12) read: 



G 2/3 (2a+m)- 1 / 3 < \a-a+\, 
G 2 / 3 (2a_m)- 1 / 3 < \a - a_|, 



(3.16) 



where the first one corresponds to a+ and the second one 
to a_. So it is a sufficient condition that 



(GA) 2 « 



4m 4 (l + m) 2 
1 (2 + 3m) 3 



(GA) « 



2 4m 4 (l — m) 2 



(-2 + 3m) 3 



(3.17) 



to conclude the existence of values of a such that there 
exists an overlapping between the WBK and the linear 
approximations. This overlapping depends only on the 
amount of radiation present in the FRW universes when 
it can be described semiclassically A <C G _1 . 

In the case of a negative cosmological constant A < 0, 
we have to deal with a unique turning point a p . Using 
similar methods to the case A > 0, we obtain that a 
sufficient condition for the existence of value of a such 
that the linear and WBK approximation hold is: 



(GA) « 



2 _ 4m 4 (m- l) 2 



(3m -2) 3 



(3.18) 



which is the case when the maximum value of the radius 
of the collapsing closed FRW universe, a p , is large enough 
(equivalently the parameter m is large). 



IV. MATCHING CONDITIONS 

Using the fact that there are values of the scale fac- 
tor, a, in which the linear and WBK approximations 
hold, we will connect the wave function on the differ- 
ent WBK regimes through the linear ones around the 
turning points, a_, a+ and a p . 



A. Positive cosmological constant 

For the case of a positive cosmological constant there 
are three WBK regimes, one corresponding to values of 



the scale factor such that a < a- 
\E'i, can be written as 



and the wave function, 



*i(a,/n) = Cexp 
+ Dexp 



+oo 

-^i(«)--^Sl(a)^ 

71=1 

1 1 + °° 



71=1 



+ 



(4.1) 



a second one, which is classically forbidden, a_ < a < o+, 
and for which the wave function, <]/„, can be expressed 
as follows 



*n(a,/„) = Aexp 
+ Bexp 



1 1 + °° 



n=l 
+ oo 



+ 



1 1 



G 



71=1 



(4.2) 



and finally, a third one in which the scale factor, a, is 
larger than a+. In this regime the wave function is 



*m(a,/n) = £exp 
+ Fexp 



i +°° 

gsrw^E?^ 

71= 1 
1 + °° 



+ 



71=1 



(4.3) 



In all these expressions, the functions with the subscript 
represent the background part of the wave function and 
satisfy the differential equation (3.10), while the func- 
tions with suffix n, n > 1, correspond to the perturba- 
tions of the wave function of the universe and they are 
solutions of the differential equation stated in expression 
(3.11). Outside the potential barrier V(a), the functions 
Sq and S™ are related to the outgoing modes and the 
functions and R™ correspond to the ingoing modes. 
On the other hand, under the potential barrier Sq and 
Rq represent growing and decreasing background terms, 
respectively, of the wave function ^ n . 

To connect the wave function of the FRW universe 
Eqs. (4.1,4.2,4.3) through the linear regimes around a_ 
and a+ we will consider that we have a unique mode /„ 
of the massive scalar field The general case can be 
easily deduced from this one. Near the turning points, 
the wave function of the closed FRW universe filled with 
radiation and a positive cosmological constant in addi- 
tion to the vacuum fluctuation of a massive conformally 
coupled scalar field, can be expressed as 

*»(a,/n) = E hmiAi(zi) + 6 mi Bi(zi)} x 

m«=0 

x exp(-y 2 t /2)H mi (y„ i ), (4.4) 
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where 



Zi 



( - iy+1 (£) 2/3 [(-^ +1 9aVo(a T )] 1/3 (a a T ) 



2G 



1/3 



+ 2 ( — ) [(-iy +1 d a Vo(a T) 



-2/3 



ft,,, 



/3„, = -(2mi + l)y/V n (a±)/2; m l e N, 
y n< = (2K(a±)) 1/4 /n- 



(4.5) 



in this expressions the index i may takes the values 1 
or 2 which correspond to the wave function on the lin- 
ear regimes around a_ and a + respectively, the potential 
Vo(a) was defined in Eq. (3.8), r ) mi and S mi are constants, 
the functions Ai(x) and Bi(x) are the Airy functions and 
H mi (x) are the Hermite's polynomials [16]. 

Let us now begin connecting the WBK wave function 
under the barrier, \I> n (a, /„), around a_ with the linear 
regime using the fact that there are values of the scale fac- 
tor a, such that both approximations, the WBK and the 
linear one, hold as was explained in Sec. IIIC. For these 
values of the scale factor, the wave function \P n (a, /„), 
using Eqs. (3.10,4.2), can be expressed near a- as 



*n(a,/„) = A exp ^VdaVo(a-) (a-a_) 3/2 



x ex P (--^(a)^) + 



B exp 



--VdaVo(a-) (a-a_) 3 / 2 



x exp(-±i£( a) /*) 



(4.6) 



Let us do some remarks about the behavior of the vari- 
ables Zi. These variables are related to the different in- 
dices rrij of the Hermite's polynomials, through the con- 
stants j3 ni Eq. (4.5). Consequently, the functions Ai(zi) 
and Bi(zi) can not be factor out from the sums over the 
different indices m^ which define the wave functions 
near the turning points. However, since we are working 
on the semiclassical framework the second term, which 
is proportional to (3 ni , on the definition of Zi is much 
smaller than the first one, as the first one is proportional 
to G~ 2 / 3 while the second one is proportional to G 1 ^ 3 . 
So the variables zi in this regime effectively do not de- 
pend on the indices rrii and the Airy functions can be 
factor out from the sum given in Eq. (4.4). Considering 
our last statement, we will connect the wave function on 
the linear regime with that of the WBK regime for values 
of zi such that the Airy functions can be approximated 
by their asymptotic behaviors (zi — -> ioo) 1 . 



Near the turning point a_, the variable z\ is positive 
under the barrier and reaches values large enough to use 
the asymptotic behavior of the Airy functions. Once the 
Ai(z\) and Bi(z\) have been substituted by their asymp- 
totic behavior in expression (4.4) for i = 1, comparing the 
resulting expression for the wave function of the universe 
*i(a,/„) near the turning points a- for a > a_ with 
the wave function *„(a, /„) in Eq. (4.6), and imposing 
the continuity of the FRW wave function, we see that 
the growing term in >]>„, related to Sq, overlaps with the 
Bi(z\) terms on the expression of \I>i, while the decreas- 
ing term in \I/ n , related to iJfj, overlaps with the Ai(zi) 
terms on the expression of Wi. Also we obtain the fol- 
lowing equations 



{ y \ 

Tmi exp ^ — |ij H mi (y ni ) 



mi— 



Bcxp 



v 2 



^2 <5 mi exp \-^y ) H mi (y ni ) = 

mi=0 



A exp 



\si{a-)fl 



(4.7) 
(4.8) 



Using the orthogonality relations of the Hermite polyno- 
mials [16] and the following formula [17]: 

/ cxp(~x 2 )H 2m (xy)dx - V^ K —r(y 2 - 1)"\ (4.9) 
we have 

= f B(l - fc_) h /{2 2l H 1 \k l l +1/2 ) 1 mi = 2h, hen 
7mi \ 0, mi = 2/i + 1, h G N, 



s _\ A{1- fc_) il /(2 2il /i!fc'! +1/2 ), mi = 2Zi, ligN 
" 0, mi = 2h + 1, he N, 



where 



2fc_ = 



y/2V n (a-) 



+ 1, 2fc_ = 



g(0 
y/2V n (a-) 



1. (4.10) 



The last expressions determine the linear behavior of the 
wave function near the point a_ in terms of the WBK 
wave function under the barrier V(a) (sec Eq. (2.3)) 
which can be explicitly seen through the dependence of 
the coefficients S mi and j mi in A and B respectively. 



To use the asymptotic behavior of the Airy functions on the 
expression (4.4), it is necessary to check that the condition \zi \ — » 
+oo, the linear and the WBK approximation overlap for some 
values of the scale factor, a, near the turning points. Indeed this 



is the case because the condition |zj| — » +oo and the validity 
of the WBK approximation (3.12) coincide near all the turning 
points a—, a_|_ and a p , in particular for the quantities of radiation 
that we are considering. 
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Now, using on the one hand the explicit expression of 
^i(a,f n ) around a_ for values of the variable Z\ such 
that z\ — > — oo and on the another hand the behavior 
of the wave function vD^ of the FRW universe outside 
the potential, V(a) (see Eq. (2.3)), ^(a, /„) near the 
turning point a_, we can obtain similar relations to the 
ones expressed in Eq. (4.8) for the coefficients j mi , 5 mi , 
C and D : 



Y (7mi/(2«) - <5 mi /2)exp(-ij)exp 



mi— 



y 2 ni 



xH mi (y ni ) = Cexp 



Y (5 mi /2 + 7 mi /(2i))exp(i-)exp( — ^ 

mi=0 



xH mi (y ni ) = Dexp 



-2<(°-)/n 



(4.11) 



where we have used the continuity of the wave function of 
the FRW universe. Finally using expressions (4.8,4.11), 
we deduce 

C - [A/2-B/(2i)]exp(-i|), 

£> = [A/2 + B/(2i)]exp(i£) , 
S^(a_) - S>_) - i£(a_) - iZ>_), (4.12) 

which determine the wave function of a FRW universe 
filled with radiation outside the potential barrier V(a), 
in terms of the wave function of a FRW universe inside 
the potential barrier V(a) and viceversa. 

Once we have obtained the matching conditions for 
the wave function in the case of a positive cosmological 
constant for the turning point a_, let us deal with the 
matching conditions for the turning point a+. For this 
purpose, we approximate the wave function, under the 
potential barrier V(a), ^ n (a, /„) defined in Eq. (4.2) near 
the point a + by 



*n(a,/„) 



3tt 
2G 



= A exp 
x(a+ -af' 2 
+-Bexp 
x(a+ -af' 2 



I + VVo(ajda - 

J a_ 

exp (-l<%(a + )f£) + 



V-d a V (a + ) x 



3tt 
2G 



^V (a)da + — y/-d a Va(a+) x 



exp --K{a + ).f 2 n 



(4.13) 



In the next step, we will match the wave function 
\& n (a, /„), using the last equation, with the wave func- 
tion on the linear regime around a +1 ^ 2 (a, f n ) expressed 
in Eq. (4.4) for i = 2. Similarly to our procedure for the 



matching conditions in a_, we use the asymptotic behav- 
ior of the Airy's functions in the expression of ^ 2 (a, fn) 
for values of the scale factor, a, such that the condition 
z 2 — > +00, the WBK approximation and the linear one 
hold and we obtain the following relations between the 
coefficients j m2 , S m2 and the constants A and B 

2 7 ™ 2 exp (-^J ff ro2 (y„ 2 ) = (4.14) 
|?7 / V^o(a)do ex 

*/ (2 — 

Y S ™ 2 CXP (^"^J H m 2 (yn 2 ) = 



m 2 =0 

A exp 

+00 



mi— 

£>cxp 



exp 



-2 R n( a +)fn 



Using as before the orthogonality relations of the Her- 
mite polynomials, we have 

f AA(1 - k + ) l -/{2 2l H 2 \k + h+1 ' 2 ) 1 m 2 = 2l 2l 
7m2 1 0, m 2 = 2l 2 + 1, 



e _ / BA-^l - fc+) i V(2 2 '^ 2 !^ +1/2 ), m 2 = 2/ 2 , 
m2 "\0, m 2 = 2? 2 + 1, 



where 
2fc+ 



V / 2K l («+) " V2V„(a+) 
3tt 
2G 



exp 



/ + VVo(a)rfa 

■/ a- 



,i 2 £N. (4.15) 



These expressions define the behavior of the wave func- 
tion v P 2 (a, /„) in the linear regime around the turning 
point a+ in term of the WBK wave function ^> u (a, /„) 
under the potential barrier V(a). 

Using the asymptotic behavior of the wave function 
^ 2 (a, /„), this time for z 2 — > +00, we can match the wave 
function of the linear regime around a+ with the wave 
function ^ m (a, f n ) outside the barrier of potential V(a). 
The continuity of the wave function implies 



+00 



Y (-7m 2 /(2i) + <W2)exp(-i^)exp(--^ 



m 2 =0 

x H m2 (y n2 ) = Eexp 

+ OC 



-^{a + )f n 



Y (^ 2 /2 + 7ro2 /(2i))exp(4)e X p(-^ 



m 2 =0 

x H m2 (y n2 ) = F exp 



(4.16) 
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and therefore 



E = 



+ 



F = 



A 



3tt 



2T XP 2G 



2i 



exp 





+ 



exp(-i-), 



exp(«^), 



S£(a+) - S™(a+) = i£( a+ ) = JC(a+)-(4-17) 

These equalities, together with Eq. (4.12) are the match- 
ing conditions for the wave function of a FRW universe 
filled with radiation, a positive cosmological constant and 
the vacuum fluctuations of a massive conformally coupled 
scalar field. Apart from these matching conditions, there 
are other conditions which ensure the good behavior of 
the wave function. These are the regularity conditions 
[8] 



RelS^a)}, Re[R 1 n (a)}>0, 
ife[S»], Re[R^(a)]>0, 
Re {S™(a)} , Re [R™(a)} > 0, 



for a < a_, 

for a_ < a < a+, 

for a + < a. (4.18) 



B. Tunneling boundary conditions of the universe 

As an example that illustrates the applications of these 
boundary condition, let us discuss the tunneling bound- 
ary conditions of the universe [5]. For these boundary 
conditions, outside and far from the potential barrier 
V(a), i.e. in the classically allowed region and for val- 
ues of the scale factor much larger than the turning point 
a + , the wave function of the universe should contain only 
outgoing modes. That is, the coefficient F (see Eq. (4.3)) 
must be equal to zero, so that no ingoing modes appear 
on the asymptotically de Sitter region. Once the tun- 
neling boundary conditions have been applied, we can 
deduce the linear combination of the growing and decay- 
ing terms that define the FRW wave function under the 
barrier, i.e. the relationship between the constants A and 
B (see Eq. (4.17)) 



\A/B\ = exp 



3tt r+ 

~Gj a _ 



(4.19) 



in agreement with the results obtained in Ref. [15] for an 
analogous system. 

So, both growing and decaying background terms are 
present in the expression of ^(a, f n ) under the barrier 
V(a). Nevertheless, we see that the growing term asso- 
ciated to 5q (a) is multiplied by the constant A which is 
exponentially smaller than the constant B, which multi- 
ply the decreasing term on ^> n (a, /„). We see that, even 



if it allows the appearance of a growing term on the clas- 
sically forbidden region, under the barrier, it is expo- 
nentially reduced. On the other hand the wave function 
^i(ffl) fn) defined on the classically allowed collapsing re- 
gion a < a_, will be a combination a of ingoing and 
outgoing modes. 



C. Negative cosmological constant 

In the case of a negative cosmological constant, A < 0, 
there is a unique turning point a p . This value of the scale 
factor separates a classically allowed region, a < a p , from 
a classically forbidden one, a p < a. The matching con- 
ditions for the wave function around a p can be deduced 
carrying out a similar analysis to the one presented pre- 
viously for A > 0. We summarize our results in what 
follows. 

In the classically allowed region (a < a p ), we will de- 
note the wave function by 



+ B A exp 



1 + °° 



n=l 



+ 



G 



n=l 



(4.20) 



while in the forbidden region (a p < a) , the wave function 
will be 



* F (a, /„) = exp 



+ oo 

7^o F («)-^E^( a )/" 



n=l 



• (4.21) 



Where we have considered only the decreasing wave func- 
tion for the asymptotic region a p < a under the barrier. 

In the linear regime, around the scale factor a = a p , the 
wave function satisfies Eq. (3.13) and can be expressed 



as 



*p(a,/n) = X] bm p Ai(z p ) +6„ lp Bi(z p )] x 



m p -0 



x exp(-yl /2)U m (y n ), 



(4.22) 



where 



3tt 



2/3 



) [d a Vo(a p )] 1/3 (a- a p ) 



H^) 1/3 [d a V(a p )}- 2 ^(3 np 



2 We consider a unique mode f„ for the massive scalar field as for 
the case of a positive cosmological constant. For the general case 
(multiples modes of the massive scalar field) the results can be 
easily generalized. 
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(3 np = -(2m p + l)^V n (a p )/2;m p e N, 

y np = (2KK)) 1/4 / n , (4.23) 

We match the wave function in the linear regime (see 
Eq. (4.22)) with the wave function under the barrier 
^F{a,f n ) (see Eq. (4.21)), taking into account that the 
background part of $ F (a, /„) is a decreasing function of 
the scale factor a, whose exponent can be approximated 
near a p by 

S F (a)=7T^d a V (a p )(a-a p )^ (4.24) 

The background part of the FRW wave function outside 
the barrier \I/ A in the neighborhood of a p corresponds to 
ingoing and outgoing modes. Therefore Sq and Rq will 
have the form 

S*£(a) = -R^a) = nr\/d a Vo(a p )(a p - a)i . (4.25) 

Taking into account that there exists values of the scale 
factor close to a p for which the linear and WBK approx- 
imation and the asymptotic condition \z p \ — > +oo hold 
simultaneously, we conclude that 

= / B A (1 - k p ) l ^/(2 2l H p lk p l - +1/2 ), m p = 2l pi l p G N, 

m p = 2l p + 1, l p e N, 



0, m p e N, 



where 



2k p — 



^2V n {a p ) 



1. 



(4.26) 



(4.27) 



Matching now the wave function *& p (a, /„) defined in the 
linear regime with \I/ A (a, /„) corresponding to the wave 
function outside the potential barrier we have 



* A (a,/„) oc exp 



i exp 



ra=l 



(4.28) 



where the proportionality symbol is related to a normal- 
ization constant that we will disregard, and the pertur- 
vative parts of the WBK wave function must satisfy 



Sn( a p) = S'n( a p) = K( a p)- 



(4.29) 



On the other hand, similar to the case of positive cosmo- 
logical constant, the functions S^(a), i?£(a) and R'n(a) 
have to satisfy regularity conditions which ensure the 
good behavior of the wave function on the WBK regime: 

Re [S*(a)] , Re [i?*(a)] > 0, for a < a p , 

Re [Sl(a)} > 0, for a p < a. (4.30) 



V. BACKGROUND WAVE FUNCTION AND 
MATTER FLUCTUATIONS 

As we saw in Sec. Ill B, the behavior of the wave func- 
tion in the WBK regime is determined by background 
and perturbative contributions (the matter associated to 
the vacuum fluctuations of a massive scalar field confor- 
mally coupled to gravity), which satisfy Eqs.(3.10,3.11). 



A. Positive cosmological constant 

In this case, there are two classically allowed regions 
and a forbidden one, when the amount of radiation 
present in the universe does not excede the maximum 
of the potential V(a). In the region a < a_, the ingo- 
ing and outgoing background parts of the wave function, 
related to Sq and Rq, can be deduced straightforwardly 
using Eq. (3.10): 

S^a) = -W = iy|" ^Voda 

ITT 1 

= y ^sa+lE^aO- mF(£ : , <*:)]- 



3m+ 1 



I a 2 (at - a 2 ) 



where 
£i = arcsin 



1/2' 



1 



OLi = 



1 + m' 



(5.1) 



(5.2) 



and F(£i, ai) and E(£ I; a r ) are the elliptic integrals of the 
first and second kind respectively [16, 17]. 

In order to study the perturbations, we will now intro- 
duce the Lorcntzian conformal time fj through 



dSo 
da 



dR'p 
da 



. 37r da 



(5.3) 



The differential equation satisfied by the perturva- 
tive parts wave function, related to and R} n , (sec 
Eq. (3.11)) can be linearized introducing the functions 
vL defined as 



Kin) = -i«y(-v)K(-fi), 



(5.4) 



where the prime denotes derivative with respect to the 
Lorentzian conformal time. In term of the functions v l n , 
the differential equation that satisfy and R\ reduces 
to 

K)"+[(n + l) 2 + M V(r?)K = 0, (5-5) 
where n is the mass of the scalar field $, and a(fj) is 



a(fj) 2 = a 2 _cd 2 



1 + m _ 1 — m 



(5.6) 
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2 -K(m),0] and cd 



is a Jacobian el- 



7 ' G \~\fl+v ""J ' ' l + m 

liptic function [16, 17]. Eq. (5.5) is a generalized Lame 
differential equation [18]. This can be seen taking into ac- 
count the explicit expression of the scale factor a(fj) given 
in Eq. (5.6), the relation cd [x + K(m), m] = — sn[x,m] 

[16], and introducing a new variable u = ^^fj, so that 

Eq. (5.5) becomes the generalized Lame equation 



d 2 < 
du 2 

with 



{N(N + l)k sn 2 [u - K(k),k] - h) v l n , (5.7) 



N(N+1) = -^-, k = \ 

A l + m 



m , h = 2<£±ll. (5.8) 



l + m 



In the classically forbidden region, (a_ < a < a+), the 
functions Rq and Sq are 



S£(a) = -i$(a) 



3tt 



V Vo da 
„, a„) - (1 - m)F(£„, a,,)] + 



(5.9) 



where 



arcsm 



1 



1/2' 



2m 



l + m 
(5.10) 

On the other hand, the differential equation satisfied 
by the perturbations S" and i?" can be simplified as be- 
fore: 



(0"-[(r* + i) 2 + A 2 WK = o, 



(5.11) 



where, now, 77 is the Euclidean conformal time related 
to the growing and decreasing background terms of the 
wave function ^ n (a, /„) by 



dSS 
da 



dR% 
da 



3ir da 



(5.12) 



a(j]) is the Euclidean scale factor given in Eq. (2.5), the 
functions z/° arc defined as 



KM 



-KY(v)K(v), (5.i3) 



and the prime denotes derivative with respect to the Eu- 
clidean conformal time 77. As the functions v l n , their ana- 
logues defined under the potential barrier V(a), sat- 
isfy also a generalized Lame equation. 

Finally, for the asymptotically de Sitter regime (a > 
a + ), the background parts of the wave function \& III (a., /„) 



are 



= 1 



+ 



a + [E(£„ I; a m ) - mF(£ m , a m )] + 
/(a 2 -a 2 )a 2 l , 

Aa -— ^- w „ 2 A h ( 5 - 14 ) 



.7T 1 

2 7x 



where 

Cm 



arcsm 



2 2 
a z - aq_ 

a 2 - a 2 . 



1/2 



m 



l + m 



(5.15) 



The perturbative parts of ^ m (a, fn) can be obtained fol- 
lowing the same procedure as for \I/i(a, /„) and ^n(a, /„). 
We first introduce the Lorcntzian time, 77, by 



dR" 1 



. 37r da 



da 



da 2 aVy' 



(5.16) 



and define the functions i+ : by 

which satisfy 

KT+[(n+l) 2 + A1 V(r7)]C = 0, 



(5.17) 



(5.18) 



where the explicit expression of the scale factor a(rj) in 
the asymptotically de Sitter regime is [17] 



a(ri) 2 = a+dc 2 



l + m _ 1 — m 



l + m 



(5.19) 



being a Ja- 
v™ also sat- 



with 77 e [0, v /i±niK(i+^)] and dc 

cobian elliptic function [16, 17]. As before, v™ 
isfics a generalized Lame differential equation. 

To obtain the explicit expression of the perturbative 
parts of the FRW universe, in the case of a positive cos- 
mological constant, it is necessary to solve the differen- 
tial equations (5.5,5.11,5.18) that satisfy the functions 
v l n , and v™. Nevertheless, it is enough to solve only 
one of them since the dependence of the scale factor a on 
the Lorentzian time fj for both classically allowed regions 
can be deduced by performing an analytical continua- 
tion of a, under the barrier V(a), from the conformal 
Euclidean 77 to the Lorentzian one 77. The differential 
equations (5.5,5.11,5.18) are related by these analytical 
continuations and so their solutions can also be related 
in the same way. Finally it must be pointed out that the 
boundary conditions that v l n , and v™ satisfy are given 
by the regularity conditions (4.18) which ensure a good 
behavior of the wave function of the universe. 



12 



B. Negative cosmological constant 

While for positive cosmological constant, A, a FRW 
universe filled with radiation can present two classically 
disconnected regions, for negative A there is just one clas- 
sically allowed region. This section is devote to the latter 
case, presenting at the end of our calculations an ex- 
plicit example in which the description of the perturba- 
tive parts of the wave function can be carried out an- 
alytically. Similarly to the preceding case, A > 0, the 
nonperturbative parts of the wave function can be ob- 
tained from expression (3.10). For the classically allowed 
region (a < a p ), the function 5* and Rq related to the 
background action are 



OA 



— — i?n 



,3tt 



y—Vo da 



1 + 3m 1 + rn 

-E(r A , s A ) + — r — F(r A , s A/ 



2A 



^-\a 2 (a 2 p -a 2 ){a 2 -a 2 n )) 



(5.20) 



while for the classically forbidden region (a p < a), the 
action can be expressed as 



So 



3tt 
Y 

IT 

2 



/ v% 



da 



' ' 5m -E(r F ,s F ) - i— — F(r F ,s F ) 



2A 

7TTO 



2A 



+ V-A— -—(a 2 -^) 



m+1 



(5.21) 



where a 2 = (1 + m)/(2A) and 



r A = arcsin 



r F = arcsin 





a 2 




la 2 - 


a 2 



a l - a l ' 



a l - a l ' 



(5.22) 



The vacuum fluctuations of the massive scaler field, $, 
yield to pcrturbative parts, S„, R„ arj d S^, in the wave 
function as described before. The analogy between the 
differential equations that govern the perturbative parts 
when containing a positive or a negative cosmological 
constant A suggests the introduction of the Lorentzian 
conformal time fj and the euclidean one rj to linearize 
Eq. (3.11) for the functions S^, R„ and S£. So, for values 
of the scale factor, a, smaller than the maximum radius 
of the collapsing FRW universe, a p , we define fj as 



dSl 
da 



da 



. 37r da 
Tdfj 1 



while for a > a p , f] is given by 



dsi 

da 



37r da 



(5.23) 



(5.24) 



Similarly to the positive A case, we introduce the func- 
tions related to and R*, in the classically allowed 
region: 



K(v) = -i{<)'{-n)/<{-nl 



(5.25) 



where the prime denotes derivative with respect to r\. In 
term of the new functions v*, Eq. (3.11) reads 



K)"+[(n+l) 2 + ^a 2 (fj)]K = 0, 



(5.26) 



and the explicit expression of the scale factor a(fj) was 
given in Eq. (2.4). For convenience, we rewrite the last 
equation in terms of the Weierstrass function V(x\w,u>') 
[16, 19] as a generalized lame equation Eq. [18] 



{<)" + 



(n + l) 2 + 



3A A 



■V(r]\iJ A ,uj' A ) 



(5.27) 

where the so called half periods, uj a and lu' a , of the Weier- 
strass function V{fj\u> A ,oj' A ) are 



w A = K 



in 



J A = iK 



rn - 



2m 



l^frn. 



Under the potential V(o), Eq. (2.3); i.e. for a p < a, 
the wave function must decreases and the unperturbed 
Euclidean space-time corresponds to an asymptotically 
adS wormhole [10]. The linearization of Eq. (3.11) for 
the functions can be made as in the preceding cases, 
that is, introducing new functions v ¥ n given by 



(5.28) 



where the prime denotes derivative with respect to the 
conformal Euclidean time and the functions v Y n satisfy 



(iO"-[(n+l) 2 + M 2 a 2 fa)K = 0. 



(5.29) 



The explicit expression of the wormhole scale factor was 
given in Eq. (2.6). In terms of the Weierstrass function, 
this equation has the form 



«)'' 



(n + I) 2 + - + ^fI^f, <4) 



0, 



3A A 

(5.30) 

where the half periods of the Weierstrass function V are 



K 



1 



2m 



/s/m, uj' f = iK 



m 



1 



2m 



/Vm. 



Summarizing, we have presented a method to deal with 
the behavior of the wave function of a FRW universe with 
a cosmological constant and filled with radiation under 
the presence of vacuum fluctuations of a massive scalar 
field conformally coupled to gravity in the semiclassical 
approximation. 
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An explicit example 



where £ F = ((uj f \lj f ,uj f ), which imply 



We will illustrate the analysis above by studying the 
stability of a radiation-filled FRW universe with nega- 
tive A in the case in which the mass of the scalar field 
is /j 2 = — 2A, with A <C lp 2 - This simple choice for the 
value of the scalar field mass allows us to solve analyti- 
cally the differential equations (5.27,5.30) for the pertur- 
bations since, in this case the generalized Lame equations 
reduce to Lame equations, whose solutions are known. 

Under the potential V (a), the perturbative parts of 
the wave function, S^, were expressed in term of the func- 
tions v^iv) m Eq. (5.28) where the functions z^(?i) satisfy 
a generalized Lame differential equation (see Eq. (5.30)). 
In the case under consideration, /i 2 = — 2A, this equation 
becomes 

«)" -\h + N{N + l)V(f] + uj f \oj f , <4)] < = 0, (5.31) 



F n = 2i [W F ((Z F \U) F ,U)' F ) - Z F £_ F ] . 



(5.35) 



As can be seen from this expression, depends on the 
parameter z F defined in Eq. (5.33). So, we have to obtain 
the possible values of z F in order to characterize F^ and 
the behavior of the linear independent solutions v\ n {r\) 
and v\ n {n). 

Since the parameter z F is defined implicitly through the 
Weierstrass function V(x\u F , uj' f ), we restrict its values to 
the fundamental rectangle [16], whose vertices coincide 



with the values 0, 



U) F , lo f , 



and to F +Wp. That is, z F can 



belong to the following ranges 

Range A : z F = f3 F , < (3 F < w F , 

Range B : z F = u> F + iS F , < S F < \ui 

Range C : z F = u' F + (3 F , < j3 F < ui F , 

Range V : z F = i5 F , < 5 F < \co' F \. 



with TV = 1 and h = (n + l) 2 - §. Since N e N this 
is a Lame equation whose solutions can be expressed as 
linear combinations of the linearly independent solutions 
vlnii) and vl n {rf) given by [18] 



a{r\ + uj f + z f \uj f , uj' f ) 



v 2n 



(v) 



a(rj + uj F \uj F ,uj F ) 
cr(—r] + uj f + z f \uj f , oj' f ) 
a(-7] + iv F \u F ,oj F ) 



exp \-i]((z f \uj F7 uj f ) 



exp [r)((z F \u> F ,i 



(5.32) 



where a{x\uj Fl u' F ) and ((x\co f ,uj f ) are Weierstrass func- 
tions [16, 19] and the parameter z F is implicitly defined 

by 



V{z f \lo f ,u' f ) = (n + lf 



(5.33) 



The differential operator that defines Eq. (5.31) is a 
Schrodinger operator whose potential is periodic as it 
can be expressed in term of the Weierstrass function 
V(x\lu f ,uj f ). Therefore, the solutions will present an 
infinite number of forbidden and allowed bands known 
as Floquet bands and the linear independent solutions 
l/ in( ? 7) and ^niv) wm be characterized by a Floquet in- 
dex F^, independent of i], such that 



"IniV + 2^f) 
V2n(V + 2u F ) 



exp(«F^K„(?7), 
exp(-iF^K„(?7). 



(5.34) 



So, for the allowed bands, defined by real values of F^, the 
amplitudes of ^(77), will be in principle bounded from 
above, while for the forbidden bands, i.e. for complex 
values of F^, the solutions will be exponentially increas- 
ing or decreasing. In the case under consideration, the 
explicit expression of the Floquet index can be deduced 
using the following propriety [19]: 

a(x + 2uj f \u! f , u' F ) = — <t(x\oj f , ut' F ) exp [2(x + w f )£f] 



When z F takes values in each of the four preceding ranges, 
its definition given by Eq. (5.33) implies [16] 



z F e Range A 
z F e Range B 
z F e Range C 
z F e Range V 



2*±i < {n+l) 2 , 
1 < (n + 1) 2 < 
^ < (n + 1) 2 < 1, 
(n + 1) 2 < 1^2. 



Now, remembering that m > 1 and n > 1, we conclude 
that the parameter z F cannot take values in the ranges C 
and T>. In the two remaining ranges A and B, the Floquet 
index can be expressed in term of Theta functions [19], 
allowing us to conclude that the range A is a forbidden 
band, i.e. F^ is complex, while the range B is an allowed 
one, i.e. F^ is real. 

Finally, the general solutions to Eq. (5.31) will be a 
linear combination of v\ n and v\ n . Using Eq. (5.28), we 
can write the functions S Y n in terms of v\ n and v\ n with 
just one free integration constant (this was expected since 
satisfies a first order differential equation). Therefore, 
for our purposes and without loss of generality we can 
write 



'In 



(5.36) 



The constants A F n have to be chosen so that the regular- 
ity conditions for hold and the functions S?^ coincide 
in the turning point a p with their counterparts in the 
classically allowed region Sl^ and R^. 

The functions ^(77), related to and R^, can be 
similarly deduced, since the differential equations (5.27) 
also reduce to Lame equations with the following linearly 
independent solutions 



a(fj + m' A + z a |cj a ,^) 
a(fj + u;' A \uj A ,L0' A ) 

a(-fj + u)' A + z A \u A ,u' A ) 
a(-f] + uj A \w A ,w A ) 



exp [-t]C(z a \oj a ,uj a )] 
exp [fj({z A \u A , uj a )] 



(5.37) 
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where the parameter z A now satisfies the following rela- 
tion 

V{z A \oj A ,J A ) = -{n+lf + ^. (5.38) 

These new functions, Vi n {fj) and V2 n (fj), are quasipcri- 
odic like their conterparts v\ n {j]) and ^2„(fy), i- e - they 
satisfy a relation analogous to Eq. (5.34) where the new 
Floquet index reads 



F n = 2i [waC(^a|w a ,w1) - ZaCa] 



(5.39) 



The values of the parameter z A can be reduced to the fun- 
damental rectangle of V(x\oj a , w' A ), but, as before, owing 
to the presence of radiation in the FRW universe (to > 1) 
and the vacuum fluctuations of the massive scalar field 
(n > 1), they can only belong to the following ranges: 



J Range £ : z A 
\ Range T : z A 



for which 



z A e Range £ 
z A e Range T 



+ Pa, 0<f3 A <u A , 
k, < S A < Ioj'A, 



1 < (n + 1) 2 < 2^1, 
< (n+1) 2 . 



Note that there is a correspondence between these in- 
equalities, satisfied in the ranges £, and T and those 
satisfied in the ranges B and A respectively under the 
barrier. This correspondence is due to the fact that, since 
the analytic prolongation of z A in Eq. (5.33) is z A = iz F 
[16], for each value of z F belonging to the range A or B, 
there is a unique value of z A belonging to the range T or 
£ respectively. 

The expression of in terms of the Theta functions 
[19] allow us to conclude that the range £ corresponds 
to a forbidden band while the range T corresponds to an 
allowed one. 

Finally, the perturbed part of the wave function, 
can be obtained using Eq. (5.25) and the solution 



(5.40) 



The function i?* can be similarly obtained provided that 
77 and A* s are substituted by — fj and A* r respectively. 

The matching conditions at the turning point a p re- 
quire that S^, S„ and i?* be equal at this point, so that 
the independent constant A F n , A^ s and A^ r satisfy the 
relation 



1 - C n + A F Jl + C n ) 
l + C n + A F n (l-C n )' 



(5.41) 



where 



c _ u; A [mg 2 (z,g F )-lngi(x, gF )]V_g_ 
ilu f [ln6 4 (x,q A ) - ln6>i(a:, g A )]' U=*a. ' 

■SWA 

q F = cxp[iirui' F / ui F ] and q A — exp[iirui' A /uj A }. 



In general, the wave function of the universe, defined 
in the whole range of the scale factor, i.e. a £ [0, +oo), 
will be a linear combination of all the wave functions 
^ a\,a 2 ... defined in Eqs. (4.20,4.21) for the allowed and 
forbidden regimes, of the form 



* = 



n 

n>l 



dAl 



U\Ai...^A\,Ai...{a,f n ), (5.43) 



where Qa\a 2 ... are the coefficient of the linear combina- 
tions and T, n denotes the set of allowed values of A F n for 
each mode, to be determined. 

The wave function 'J must satisfy a regularity condi- 
tion that in terms of the functions S^, and i?* become 
a set of inequalities requiring that the real part of these 
functions be positive as stated in Eq. (4.30). These in- 
equalities select the allowed values of A F n for each mode. 
In others words, the allowed values of A F n will be those 
for which the minimum, M[A F n ], of Re[S F n ], Re[S*] and 
Re[Rn] for all times (fj and 77) is strictly positive. Nat- 
urally, this minimum will depend both on the amount 
of radiation present in the universe and the cosmological 
constant, which are jointly represented by to, and on the 
mode n itself. This dependence is encoded in the param- 
eter z F (and its analytical continuation z A ) which may 
belong either to the range A or B. 

We have plotted the contours corresponding to con- 
stant values of the minimum M[A^] in the complex A F n - 
plane for the two possible ranges A and B as well as for 
the value z F — lo f that defines the border between both 
ranges. These contour plots are shown in Fig. 2. It can 
be seen that for z F € range A, the minimum M is positive 
outside the circle of unity radius and centered at the ori- 
gin, i.e. for \A F n \ > 1. Therefore, for fixed amount of radi- 
ation (to = 241 in Fig. 2), the allowed values of A F n with 
n such that z F € range A are £„ = {A F n e C, \ A F n \ > 1}. 
Fig. 2- (a) shows the contour plot of M for a z F in the 
range A corresponding to n = 15 and to = 241. The 
modes n such that z F e range B exhibit a more compli- 
cated behavior. Depending on the specific mode under 
study the allowed values A F n are either the upper or the 
lower complex planes, i.e. £„ = {A F n e C, Tro^J > 0} 
for some n and X„ = {A F n £ C, iro[A^] < 0} for the oth- 
ers. Fig. 2-(b) shows the contour plot of M for a z F in 
the range B corresponding to n = 5 and to = 241. Other 
modes will present either a similar plot or the mirror im- 
age with respect to the real axis, as already discussed. 
Finally, we see that for those amounts of radiation for 
which to = 2(ub + l) 2 — 1 for some ub € N, the mode 
riB such that z F belongs to both ranges A and B, the 
allowed values of A F n are £„ = e C, \ A F n \ > 1} as in 
Fig. 2-(c). Note however that unlike Fig. 2- (a), the values 
l-^nl — 1 correspond to negative infinite minima M(A F n ). 
In general, for a given amount of radiation, there will not 
exist any such boundary mode and only for very specific 
fine-tuned amounts of radiation this will happen. 

The existence of a set £„ for fixed values of the pa- 
rameter to is similar to the case studied in [7]. There 
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FIG. 2: This figure shows the contour plots corresponding to constant values of M[j4^] in the complex ^4^-plane where Af[.<4£] 
is defined as the minimum of Re[S^], Re[S%] and Re[R*] for all times. In these contour plots, the scale of grey represents 
the magnitude of M: the darkest (black) one corresponds to negative infinite values of M while the lightest ones correspond 
to positive values and therefore determine the allowed A^'s. The parameter m related to the amount of radiation and the 
cosmological constant has been set equal to m = 241. As discussed in main text, by the regularity conditions Eq. (4.30), the 
allowed values of A F n are those for which M is positive. The contour plot in Fig. 2(a) corresponds to the mode n = 15. This 
mode for the chosen value of m is such that z F belongs to the range A. The set {^15 £ C, 1 < ^15 1} corresponds to the allowed 
values of A F 1B by the regularity conditions. Fig. 2(b) represents the contour plot for the mode n — 5 for which z F belongs to 
the range B. The values of A Y 5 for which the regularity conditions hold are {A ¥ 5 £ C, 1 < ^H}. Finally, in Fig. 2(c) shows 
the contour plot of M for n = 10. In this case z F lies in the boundary of the range A and B. The allowed values of A ¥ w are 
{A\q £ C, 1 < |j4.io |} an d M reaches infinite negative values for any other value of A\ . 



the author constructed the wave function of the gravi- 
tons in a de Sitter background for different boundary 
conditions. When the decreasing wave function for the 
gravitons under the potential barrier V(a) for K = was 
picked up (boundary conditions similar to the one consid- 
ered in [9]), the wave function was not uniquely defined 
or equivalently it can be constructed as a superposition 
analogue to Eq. (5.43). The situation is rather different 
when the increasing wave function of the gravitons under 
the potential barrier V(a) is chosen: in this case there is 
a unique wave function. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the quantum behavior 
of a radiation-filled FRW universe with a cosmological 
constant in the presence of vacuum fluctuations repre- 
sented by a massive scalar field conformally coupled to 
gravity. 

In the semiclassical approximation, the wave function 
of the universe can be expressed as linear combinations 
of outgoing and ingoing modes in the classically allowed 
regions and as increasing and decreasing modes in the 
classically forbidden ones. For negative cosmological con- 
stant, the matching conditions have been deduced for 
natural boundary conditions which pick up the decreas- 
ing wave function in the forbidden region (i.e., in the 
asymptotically anti de Sitter Euclidean wormhole). For 
positive cosmological constant, the matching conditions 
have been worked out for arbitrary boundary conditions 



and have been applied to the specific case of the tunnel- 
ing boundary conditions of the universe [5] . In this case, 
the wave function describes a de Sitter-like universe that 
contains only outgoing modes in the asymptotically de 
Sitter region. These boundary conditions allow the pres- 
ence of decreasing and increasing modes under the poten- 
tial barrier. However, the ratio between the coefficients 
of the increasing and decreasing modes is exponentially 
suppressed. 

Especially important are the regularity conditions that 
we have imposed on the wave functions, namely, that 
they must be finite and well-behaved everywhere and for 
every field configurations. This conditions impose im- 
portant restriction on the allowed wave functions as we 
have shown. In particular, they guarantee that there are 
no divergences that could be interpreted as leading to 
instabilities of the background configuration (asymptot- 
ically anti de Sitter Euclidean wormhole or asymptoti- 
cally de Sitter Lorentzian region, depending on the value 
of the cosmological constant). Therefore, we have seen 
that such regularity conditions are not empty, at least 
for some values of the cosmological and the scalar field 
mass. Furthermore, we have also shown, in this case, 
that they are not too restrictive either. Indeed, there 
exist a whole sets of wave functions, characterized by a 
continuous index for each mode, which arc regular and 
therefore feasible candidates for quantum states. In this 
sense, it is worth noting that this is true not only in gen- 
eral but also for each mode separately, i.e, the regularity 
conditions allow contributions to the wave function from 
every single mode without exception. 
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With these ingredients, we have obtained explicit so- 
lutions for the background wave function and nonlin- 
ear differential equations that govern the behavior of the 
vacuum fluctuations. Appropriate linearization of these 
equations gives rise to generalized Lame differential equa- 
tions. As an application of the general procedures de- 
scribed in this work, we have fully solved the problem of 
obtaining the wave function of an asymtotically anti do 
Sitter wormhole and its quantum stability against vac- 
uum fluctuations represented by conformally coupled 
scalar field whose mass is given by \i = — 2A. This specific 
choice has allowed us to solve the generalized Lame equa- 
tions and thus fully study the quantum behavior of the 
vacuum fluctuations. As we have already discussed, the 
wormhole boundary conditions and the regularity con- 
dition that the wave function be finite for all possible 
values of the scale factor and field configurations provide 
the set of allowed quantum wormhole states, which are 
therefore stable under vacuum fluctuations. It is worth 
noting that the boundary and regularity conditions do 
not select a single quantum state as happened in Ref. [8] , 
but a set of allowed quantum states labeled by a continu- 
ous parameter for each mode. This situation is analogous 
to one of the cases studied in Ref. [7], where the autor 
obtained the wave function of the gravitons in de Sit- 
ter space. If this wave function contains only decreasing 



modes, the regularity conditions do not select a unique 
quantum state for each mode, in opposition to the case 
where the boundary condition picked up the increasing 
wave function. 

In the models considered in Ref. [7, 8], as well as the 
one studied in this paper the wave functions is well- 
behaved in the classically forbidden region, due to quan- 
tum gravity effects, indeed it is not divergent. All these 
examples show that in the classically forbidden regions, 
due to quantum gravity effects, do not led in principle 
to infiniteness of the wave function of the universe and 
that it is well-behaved in opposition to the case studied 
in Ref. [9], where catastrophic particle creation led to 
divergence of the wave function. 
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